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ABSTRACT 

Context. The Yarkovsky effect is a weak non-gravitational force leading to a small variation of the semi-major axis of an asteroid. Us¬ 
ing radar measurements and astrometric observations, it is possible to measure a drift in semi-major axis through orbit determination. 
Aims. This paper aims to detect a reliable drift in semi-major axis of near-Earth asteroids (NEAs) from ground-based observations 
and to investigate the impact of precovery observations and the future Gaia catalogue in the detection of a secular drift in semi-major 
axis. 

Methods. We have developed a precise dynamical model of an asteroid’s motion taking the Yarkovsky acceleration into account and 
allowing the fitting of the drift in semi-major axis. Using statistical methods, we investigate the quality and the robustness of the 
detection. 

Results. By filtering spurious detections with an estimated maximum drift depending on the asteroid’s size, we found 46 NEAs 
with a reliable drift in semi-major axis in good agreement with the previous studies. The measure of the drift leads to a better orbit 
determination and constrains some physical parameters of these objects. Our results are in good agreement with the l/D dependence 
of the drift and with the expected ratio of prograde and retrograde NEAs. We show that the uncertainty of the drift mainly depends on 
the length of orbital arc and in this way we highlight the importance of the precovery observations and data mining in the detection of 
consistent drift. Finally, we discuss the impact of Gaia catalogue in the determination of drift in semi-major axis. 

Key words. Astrometry - Celestial mechanics - Minor planets, asteroids: general 


1. Introduction 


Highlighted in the early 20 th century, the Yarkovsky effect is a 
non-gravitational force due to an anisotropic emission of ther¬ 
mal radiation. It mainly depends on the size, the density, the 
spin, and the thermal characteristics of the asteroid. The effect 
has a diurnal component related to the asteroid rotation, which 
is generally dominant, and a seasonal component related to the 
orbital motion ( [Vokrouhlicky et al.|2000) i. 

The Yarkovsky effect leads to an acceleration on the asteroid 
that modifies its semi-major axis as a secular drift. It has im¬ 
portant consequences for the long-term evolution of asteroids. 
In particular, the Yarkovsky effect could explain the transporta¬ 
tion of asteroids from the main belt region to the inner solar sys¬ 
tem ( Bottke et al.[2006) . It also becomes necessary for accurate 
orbit determination or prediction of trajectory of near-Earth as¬ 
teroids (NEA), in particular, in the context of close approach or 
impact predictions (Farnocchia et al.|2013a Farnocchia & Ches- 


|ley|2014[|Chesley et al.|20l4]T 

The first detection of Yarkovsky acceleration has been per¬ 
formed thanks to radar measurements of the asteroid (6489) 
Golevka ( |Chesley et al.|2003] ). |Vokrouhlicky et al.| < [2008| > mea¬ 
sured the Yarkovsky effect using only astrometric observations 
of the asteroid (15 2563) 1992BF _ 

More recently, Nugent et al. ( 201 2b|> detecte d a measurable 
secular drift for 54 NEAs and |Famocchia et al. (2013b! for 21 
NEAs using more restrictive criteria. Nugent et al. (2012b I used 
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Farnocchia et al. 


the Orbfit packagsHfor orbit determination and _ 

P013h| ) used JPL’s Comet and Asteroid Orbit Determination 
Package and cross-checked their results with Orbfit. 

In this paper, we perform an independent study of the detec¬ 
tion of the drift in semi-major axis of NEAs. We develop our 
own dynamical model allowing the orbit determination and the 
measurement of a secular drift in semi-major axis (Sect. [2}. We 
aim to detect a secular drift for the 10953 NEA0(1629 num¬ 
bered and 9324 unnumbered) from ground-based observations. 
We first verify the robustness and quality of the detection with 
statistical methods (Sects. EH & |3.2| > and finally, by introduc¬ 
ing filters to avoid unreliable detections, we detect a consistent 
measurement for 46 NEAs (Sect. [3]». As a long period of ob¬ 
servations appears necessary for an accurate detection, we anal¬ 
yse how data mining through precovery observations can help to 
refine the measurement of a secular drift (Sect. |4ji. Finally, we 
highlight the expected improvement with the Gaia mission for 
the detection of Yarkovsky acceleration (Sect. [5}. 


2. Orbit determination 

To perform the orbit determination and to detect the drift, we 
build a dynamical model of an asteroid’s motion taking the main 
perturbations into account (Sect. and we add the Yarkovsky 
effect as a transverse force depending on the orbital elements 

1 http://adams.dm.unipi.it/orbfit/ 

2 As of 12 May 2014. 
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and on a drift in semi-major axis (Sect. [2T2|> . The observations 
and their treatment are presented in Sect. [23] Finally, the fitting 
process, used to determine the state-vector of the asteroid (i.e. 
the position and the velocity of the asteroid at a specific date) 
and the value of the drift a, is presented in Sect.|2.4| 


2.1. Dynamical model 

To measure a drift in semi-major axis of asteroids, we developed 
a numerical integration of the motion of an asteroid (NIMA, 
thereafter). The equations of the motion are numerically inte¬ 
grated with a Gauss-Radau method (Everhart 1985[ >. The dynam¬ 
ical model includes: 


gravitational perturbation. The effect is considered as a trans¬ 
verse force: 


a Y = A 2 g(r) t 


(1) 


where t is the normalised transverse vector (r x v) X r with r 
and v the heliocentric position and the velocity of the asteroid, 
respectively. 

The expression g(r) is an r-function developed in |Marsden| 
|et al.|(|1973[) and A 2 is a parameter associated with this force. As 


Farnocchia et al. ( |20l 3b i, we used a simple expression with 
g(r) — ( ro/rY where ro — 1 au is a normalising parameter. 

The Gauss equations give the variation of the semi-major 
axis related to the transverse acceleration 


- the perturbations of the planets, Pluto and the Moon, whose 
positions are given by JPL Ephemeris DE421 ([Folkner et al.| 

12509 ] ); 

- the perturbations of the ten main asteroids whose posi¬ 
tions have been preliminary computed by NIMA and turned 
into Chebychev ephemeris to make quick computation. The 
masses of the main asteroids are provided in Table [T] 

- the corrections of relativistic effects of Sun, Jupiter, and the 
Earth as described in |Moyer| ( |2000) >; 

- the Yarkovsky effect (see Sect.|2.2[). 


Table 1 . Masses of the ten biggest asteroids from 
used in NIMA. 


Folkner et al. 


(2009) 


asteroid 

mass (xlO lI M Q ) 

(1) Ceres 

468.517 

(4) Vesta 

132.844 

(2) Pallas 

100.985 

(10) Hygiea 

40.418 

(704) Interamnia 

18.566 

(16) Psyche 

16.826 

(15) Eunomia 

12.342 

(511) Davida 

12.342 

(3)Juno 

11.574 

(52) Europa 

10.203 


da 2 a V 1 - e 2 A 2 r 2 

— - a - -;- 

dt nr 3 

where a is the semi-major axis, e the eccentricity, n the mean 
motion, and r is the heliocentric distance of the asteroid. 

Consequently, by averaging the previous relation on an or¬ 
bital period, the mean variation of a is 

2A 2 r 2 
na 2 ( 1 - e 2 ) 

Finally, the Yarkovsky acceleration is written as a transverse 
force only depending on the orbital elements and the drift in 
semi-major axis (now written as a): 


n a 2 { 1 -e 2 ) 

a Y = --~- (a) t. 

2 r L 


( 2 ) 


This expression is equivalent to those given in |Chesley et al.| 
(2008J and Mouret & Mignard (2011), and is consistent with 
models in Far nocchia et al.| ( |2013b| ) where they estimated the 
parameter A 2 . We can estimate either A 2 parameter or the mean 
drift in semi -major axis a. Typical valu es for the drift are 10 -3 to 


10 4 au/My (Vokrouhlicky et al. 2000: ), corresponding to a range 
of 10" 13 to 10 13 au/d for A 2 . Thereafter we use 10 4 au/My as 
a unit for the drift in semi-major axis a and its uncertainty <x„, 
and 10 13 au/d 2 for A 2 parameter and its uncertainty. 


All solar system objects in the dynamical model are con¬ 
sidered point masses. The equations of variation are also nu¬ 
merically integrated with the equations of motion to fit the 
model to observations by a Levenberg-Marquardt algorithm (see 
Sect. |2.4| ). To be consistent with the planetary ephemeris, we use 
the masses of the asteroids provided by Folkner et al. ( |2009| > 
and given in Table [I] The masses of some big asteroids are not 


Interamnia is 4.751km 3 /s 2 in Farnocchia et al. 

2013b) whereas 

the mass is 2.464km 3 /s 2 in 

Folkner et al. ( 

200S 

). Moreover, we 


compared our results for the 46 selected NEAs (see Sect. [3)4 1 


by including the 25 biggest asteroids from Folkner et al.|(|2009 


as well as the 10 biggest asteroids, and we noticed only minor 
differences. Consequently, and also for time computation, we 
choose the ten biggest asteroids corresponding to a mass greater 
than 10 _11 M o . 


2.2. Yarkovsky acceleration modelling 

We consider the Yarkovsky effect a comet-like model as pre¬ 
sented in jMarsden et al.| ( |1973| ), and we use the effect for non- 
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2.3. Observations and treatement 


The asteroid observations come from the Minor Planet Center 
databas^Jas of 12 May 2014. Observations are also available on 
the AstDyS databasq 4 ] Observations suffer from several sources 
of systematic errors that can be crucial for asteroid orbit deter¬ 
mination (Carpino et al. 2003). Zonal errors of stellar catalogues 
are one of the systematic errors for observations. Chesley et al. 
( [2010[ > proposed a method to remove the bias for several stel¬ 
lar catalogues. Bias corrections are given for each observation in 
the AstDyS database. Finally, the outliers are eliminated thanks 
to the procedure described in Carpino et al. ( j2003[ ). 


2.4. Fitting and weighting processes 

The fitting process consists in the determination of seven 
parameters (the six components of the position and velocity at a 
given date and the drift in semi-major axis d or the A 2 parameter 
defined in Sect. \2.2) that minimises the residuals (observed 
positions minus computed positions). This determination uses 

3 http://www.minorplanetcenter.net 

4 http://hamilton.dm.unipi.it/astdys/ 
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a Levenberg-Marquardt algorithm by iterative corrections of 
each parameter. For each iteration, the corrections to apply are 
determined thanks to the equations of variation and the least 
square method (LSM; for more information, see for example 
Desmars et al. 2009j l. In the LSM, a weighting matrix V 0 b s is 
required and usually considered as a diagonal matrix where the 
diagonal components are ef = 1 /err where err is the estimated 
variance of the observation i. 


In orbit determination, the main difficulty is to give an ap¬ 
propriate weight to each observation. Carpino et al. ( |2003 ! and 
Chesley et al.| ( 20T0] > discussed weighting schemes and showed 
the orbit determination improvement by weighting observations, 
according to the observatory and the stellar catalogue used for 
the reduction. 

The AstDyS database provides the weight for each observa¬ 
tion. As a general rule, we use the value given in AstDyS, except 
for specific cases. As in Farnocchia et al. ( [2013b i, we consider 
that observations before 1950 have a weight of 10 arsec, and 
observations from 1950 to 1990 a weight of 3 arcsec. These ob¬ 
servations are mainly photographic plates affected by low sensi¬ 
tivity and reduced with less accurate stellar catalogue. The other 
case is when observations have been carefully reduced again, as 
in |Vokrouhlicky et al.|(|2008|), with four precovery observations 
of (152563) 1992BF in January 1953. 


3. Detection of the drift 


In this section, we try to detect a drift in the semi-major axis of 
NEAs. We first study the consistency of the drift and its uncer¬ 
tainty determined by LSM by comparing with statistical methods 
for three representative NEAs (Sect. ED- We also look after the 
robustness of the drift we determined (Sect. [3/2) a nd we present 
the drift in semi-major axis for NEAs (Sect. 3.3| & |3.4| i and a 
comparison with the previous studies (Sect. |3T5| >. Finally, we dis¬ 
cuss the constrains in some physical parameters of the asteroids 
given by the drift (Sect. [3T6|». 


[et al.| ( 200T) and consists of adding a random noise on each ob¬ 
servation and then fitting the dynamical model to the new set of 
observations, giving a new drift value. The bootstrap re-sampling 
(BR) was developed by |Efron| ( |1979) > and consists in creating a 
new set of observations by sampling the original set of obser¬ 
vations with a draw with replacement. The new set of observa¬ 
tions has the same number of elements as the original set. Con¬ 
sequently, some observations can appear several times whereas 
others do not appear at all. The dynamical model is then fitted to 
the new set, giving a new drift value. 

For both methods, the process is repeated many times and a 
measure of the uncertainty is given by the standard deviation of 
all the drift values. We consider three significative NEAs: 


1. (2340) Hathor for which we determined a = (-14.312 + 
3.297) x 10~ 4 au/My, represents NEAs whose signal to noise 
ratio (S/N) is greater than 3; 

2. (138852) 2000WN10 with a = (15.025 ± 7.159) x 

10~ 4 au/My, represents NEAs with an average S/N; 

3. (153792) 2001VH75 with a = (12.812 ± 165.373) x 
10~ 4 au/My, represents NEAs with very bad detection of the 
drift. 

For both methods (MCO and BR), we generate 1000 sam¬ 
ples. Table [2] presents the value of the drift and the standard 
deviation obtained with the fitting process (LSM) and the two 
different methods (MCO and BR). 

The estimations of the uncertainty of the drift are equivalent 
for LSM and MCO. With the bootstrap method, we have slight 
differences in particular for NEAs with a bad detection of the 
drift. Actually, LSM and MCO make use of equivalent assump¬ 
tions (independence of observations, Gaussian errors of obser¬ 
vations) whereas the BR makes use of minimal assumptions (in¬ 
dependence of observations only). For NEAs with a good S/N, 
which is one of our criterions of selection (see Sect. |3.4| ) the 
three methods give similar results. The standard deviation <t„ de¬ 
termined with LSM is a good estimator of the drift in semi-major 
axis uncertainty and can be used as a measure of the uncertainty. 


3.1. Estimation of uncertainty 


3.2. Solution stability 


The fitting process gives the orbital elements (more specifically 
the state vector) of the asteroid and the drift in the semi-major 
axis. The LSM also provides the covariance matrix A and the 
standard deviation of each parameter is given by the root square 
of diagonal elements. Thus, for the drift in semi-major axis, the 
standard deviation is given by y/Arj. This can be considered as 
a measure of the uncertainty of the drift value cr^. According 
to LSM theory, the uncertainty of the drift only depends on the 
date of observations and the weights given to these observations. 
In particular, it does not depend on the observation coordinates 
whereas the drift value depends on the date, the weights, and 
the coordinates of the observations. This aspect is important be¬ 
cause we can study the evolution of the drift uncertainty related 
to the precision (or the weight) of the observations without the 
knowledge of the coordinates themselves and thus we can de¬ 
termine the drift uncertainty with simulated observations (see 

sect.g&|g. 

To verify if this measure is a good estimation of the uncer¬ 
tainty, we compare cy, with accuracies obtained by two statis¬ 
tical methods. In this context, we apply a Monte Carlo process 
on the observations and a bootstrap re-sampling as described in 
[Desmars et al.| ( |2009| >. The Monte Carlo method on the obser¬ 
vations (MCO) comes from a technique developed by |Virtanen| 


The stability of the solution can be tested with a random hold¬ 
out method. It consists of randomly extracting a percentage of 
the observations, fitting the model with these observations, and 
repeating the process several times. The changes in the fitted pa¬ 
rameter gives information on the stability. In this section, we use 


the same representative NEAs as in Sect. 3.1 and we randomly 
remove 10% of the initial set of observations. The process is 
repeated 1000 times. Table 3]presents the mean (/./) and the stan¬ 
dard deviation (cr) of the drift in semi-major axis a and its un¬ 
certainty cr a , obtained with random holdout method for (2340) 
Hathor, (138852) 2000WN10, and (153792) 2001VH75. 


Table 3. Statistics of drift in semi-major axis a and its standard de¬ 
viation <j ti (in 1(T 4 au/My) obtained with random holdout method for 
(2340) Hathor, (138852) 2000WN10, and (153792) 2001VH75. 


NEA 

a 

cr 

ma 

E 

cr 

(2340) Hathor 

-14.358 

0.654 

3.495 

0.129 

(138852) 2000WN 10 

15.089 

1.254 

7.473 

0.084 

(153792) 2001VH75 

14.223 

37.476 

176.367 

6.120 
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Table 2. Drift in semi-major axis and its standard deviation (in 10 4 au/My) obtained with three methods for (2340) Hathor, (138852) 2000WN10, 
and (153792) 2001VH75. 


Asteroid 

(2340) Hathor 

(138852) 2000WN10 

(153792) 2001VH75 

LSM 

-14.312 + 3.297 

15.025 + 7.159 

12.812+ 165.373 

MCO 

-14.417 + 3.294 

16.612 + 6.622 

18.183 + 158.107 

BR 

-14.497+ 1.798 

14.998 ± 3.747 

25.342+ 185.706 


The means for a and cr a have to be compared with values 
given in Table[2]with LSM. By removing 10% of the initial set, 
the values obtained for the drift and its uncertainty do not vary 
too much. It shows that the drift in semi-major axis a we detected 
and its uncertainty cr A are robust and stable, even for inaccurate 
determination of a. 


3.3. Drift determination for numbered NEAs 

In this section, we try to determine a drift in semi-major axis for 
the 1629 numbered NEAs as they have been observed over sev¬ 
eral oppositions. Figure [T| represents the length of the observed 
arc, the uncertainty of the drift in semi-major axis, and the num¬ 
ber of observations for the 1629 numbered NEAs. 



1 10 100 


observed arc (years) 

Fig. 1. Observed arc, uncertainty of the drift in semi-major axis (in 
au/My), and absolute magnitude H for the 1629 numbered NEAs. NEAs 
with radar measurements are represented with crosses and NEAs with¬ 
out radar measurements with bullets. 


For most NEAs, the uncertainty of the drift we detected is 
greater than 10 5 au/My, leading to a bad S/N and inconsis¬ 
tent detections. As the order of magnitude fo r the drift is about 
ICE 3 - ICE 4 au/My (Vokrouhlicky et al. 2000), an uncertainty of 
the same order is at least required for a consistent determination 
of the drift. 


Figure[T]shows that a long period of observations (at least 20 
years) is necessary to have a small uncertainty of the drift. Nev¬ 
ertheless, with radar measurements, the period of observations 
can be smaller because they lead to important constraints in the 
orbit determination (see for example Desmars et al.| |2013j l and 
in particular for the detection of Yarkovsky acceleration (see the 
case of (6489) Golevka in |Chesley et al.|2003[ l. As the length of 
orbital arc is longer for bright NEAs, we do not notice a clear 
correlation between H magnitude and drift uncertainty. 
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3.4. Drift determination for selected NEAs 


We try to detect a drift in semi-major axis for the 10953 NEAs 
(1629 numbered and 9324 unnumbered) as of 12 May 2014. To 
select reliable detections of drift in semi-major axis, we adopt 
two filters. 

The first filter is the signal to noise ratio defined as S/N= 
d|/(X n . jFarnocchia et al.| ( [2013b] ) used S/N> 3 as a threshold to 
define reliable detections but they also considered lower values 
for some specific cases. This threshold appears restrictive be¬ 
cause a detection with S/N~ 2 also leads to consistent drifts. In 


contrast, Nugent et al. (J2012b) used S/N> 1 ,which seems to be 
a nonrestrictive threshold. If we assume a Gaussian distribution 
for the drift value and if S/N=l, the probability of detecting a 
drift value with an opposite sign from the real value (for exam¬ 
ple a positive drift where the real drift is negative or, conversely, 
a negative drift where the real drift is positive) is about 15.9%, 
whereas if S/N=2, this probability is 2.3%. In that context, we 
adopt S/N> 2 as an appropriate filter. 

The second filter is to avoid spurious detections: for example, 
the detection of a strong drift for large NEAs. We estimated a 
maximum value for the drift in the semi-major axis according to 
the absolute magnitude H of the NEA, 


1() 0.2H-6.85 au/My 


(3) 


The maximum value assumes average values of semi-major 
axis, albedo, and bulk density for NEA according to |Chesley| 
|et al.| (2002). The full details of the estimation are given in Ap¬ 
pendix^ This parameter remains an approximation and can be 
underestimated by a factor of 3 for NEA with low density and 
by a factor of 5 for NEAs with high eccentricity ( |Vokrouhl icky 
et al. 2000[). The ratio / ; = \d\/d max can help to filter spurious 


detection. The parameter if is expected to be smaller than 1 but 
by considering the approximate value of the maximum drift, we 
select NEAs for which £ < 3. 

Finally, Table[4]presents the selected NEAs, for which S/N> 
2 and £ < 3, with the drift in semi-major axis a, its uncer¬ 
tainty <Ta, the At parameter and its uncertainty cr Al , the signal to 
noise ratio S/N, the absolute magnitude H, the diameter D (from 
EARN databas^J, the semi-major axis (in au), the eccentricity, 
the £ ratio, and the time span of observations. An additional as¬ 
terisk indicates that radar measurements are available for this 
NEA. We also add three NEAs with <Ta < 10 4 au/My because 
as discussed later, this sort of a precision indicates a drift close 
to zero, which provides an initial constraint. 

Most of £ values are close to 1, indicating that a drift de¬ 
tection is currently possible when the drift rate is close to its 
maximum value. Moreover, £ close to 2 or 3 are associated with 
a high eccentricity. 


5 http://earn.dlr.de/nea/ 
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Table 4. Drift in semi-major axis and its uncertainty, A 2 parameter and its uncertainty, S/N, H magnitude, diameter, semi-major axis, eccentricity, 
£ ratio, and orbital arc for 43 selected NEAs with S/N> 2 and £ < 3 and 3 NEAs with <r„ < 10 -4 au/My (see text). 



NEA 

a 

O'a 

^2 

O'At 

S/N 

H 

D 

a 

e 

? 

time-span 



(10~ 4 au/My) 

(10- 15 

au/d 2 ) 



(m) 

(au) 




101955 

Bennu 

-18.611 

0.093 

-44.581 

0.222 

200.6 

20.2 

484 

1.13 

0.20 

0.99 

1999-2013* 

152563 

1992BF 

-11.297 

0.839 

-23.482 

1.744 

13.5 

19.8 

510 

0.91 

0.27 

0.94 

1953-2011 


2005ES70 

-69.056 

7.823 

-120.810 

13.687 

8.8 

23.8 

- 

0.76 

0.39 

0.89 

2005-2013 

3908 

Nyx 

11.545 

1.754 

29.781 

4.525 

6.6 

17.3 

1000 

1.93 

0.46 

2.86 

1980-2014* 

4179 

Toutatis 

-2.087 

0.363 

-4.718 

0.821 

5.7 

15.3 

2800 

2.53 

0.63 

1.37 

1976-2014* 


2006CT 

-49.077 

8.607 

-114.620 

20.101 

5.7 

22.2 

- 

1.10 

0.23 

1.22 

1991-2014* 

152664 

1998FW4 

15.537 

2.910 

27.818 

5.210 

5.3 

19.7 

- 

2.52 

0.72 

1.28 

1994-2013* 

2062 

Aten 

-5.572 

1.084 

-12.469 

2.426 

5.1 

17.1 

1300 

0.97 

0.18 

1.56 

1955-2014* 

1862 

Apollo 

-2.314 

0.491 

-4.537 

0.963 

4.7 

16.4 

1400 

1.47 

0.56 

1.02 

1930-2014* 


2009BD 

-418.170 

88.640 

-986.640 

209.140 

4.7 

28.1 

4 

1.01 

0.04 

0.67 

2009-2011 


2007TF68 

-98.134 

21.758 

-254.770 

56.487 

4.5 

22.7 

- 

1.41 

0.26 

2.02 

2002-2012 

29075 

1950DA 

-2.729 

0.656 

-6.217 

1.495 

4.2 

17.6 

2000 

1.70 

0.51 

0.75 

1950-2014* 


2004KH17 

-42.423 

10.228 

-63.288 

15.258 

4.1 

21.9 

197 

0.71 

0.50 

1.26 

2004-2013* 

2340 

Flathor 

-14.131 

3.529 

-24.377 

6.087 

4.0 

20.2 

300 

0.84 

0.45 

1.00 

1976-2012 

37655 

Illapa 

-11.351 

2.925 

-14.090 

3.631 

3.9 

18.0 

- 

1.48 

0.75 

2.23 

1994-2013* 

54509 

YORP 

-36.240 

10.581 

-80.829 

23.600 

3.4 

22.6 

100 

1.00 

0.23 

0.79 

2000-2005* 

2100 

Ra-Shalom 

-4.502 

1.318 

-7.827 

2.291 

3.4 

16.1 

2240 

0.83 

0.44 

1.89 

1975-2013* 

4581 

Asclepius 

-19.683 

6.031 

-40.896 

12.532 

3.3 

20.7 

- 

1.02 

0.36 

1.00 

1989-2013* 

1620 

Geographos 

-2.708 

0.913 

-6.315 

2.129 

3.0 

16.5 

5000 1 

1.25 

0.34 

1.77 

1951-2013* 

350462 

1998KG3 

-24.853 

8.343 

-62.171 

20.871 

3.0 

22.1 

- 

1.16 

0.12 

0.68 

1998-2013 

6489 

Golevka 

-5.468 

1.796 

-13.120 

4.310 

3.0 

19.1 

350 2 

2.52 

0.60 

0.62 

1991-2011* 

283457 

2001MQ3 

-14.692 

5.027 

-40.926 

14.005 

2.9 

19.1 

500 

2.23 

0.46 

1.73 

1951-2011 

256004 

2006UP 

-95.050 

35.471 

-256.300 

95.646 

2.7 

23.0 

- 

1.59 

0.30 

1.64 

2002-2012 

363599 

2004FG11 

-28.744 

10.805 

-40.589 

15.257 

2.7 

21.0 

152 

1.59 

0.72 

1.28 

2004-2014* 

85770 

1998UP1 

-16.118 

5.982 

-33.419 

12.403 

2.7 

20.5 

- 

1.00 

0.34 

0.95 

1990-2013 


1999SK10 

-20.989 

7.650 

-52.845 

19.260 

2.7 

19.7 

400 

1.76 

0.44 

1.73 

1999-2014 


1995FJ 

42.330 

16.182 

96.248 

36.795 

2.6 

20.5 

300 

1.08 

0.27 

2.69 

1995-2013 


2001QC34 

-32.160 

12.255 

-77.598 

29.570 

2.6 

20.1 

- 

1.13 

0.19 

2.22 

2001-2014 

99907 

1989VA 

11.931 

4.520 

15.498 

5.871 

2.6 

18.0 

550 

0.73 

0.59 

2.33 

1989-2012 

164207 

2004GU9 

-67.288 

26.497 

-155.570 

61.261 

2.5 

21.1 

163 

1.00 

0.14 

2.81 

2001-2014 


1999MN 

28.367 

11.432 

30.576 

12.322 

2.5 

20.6 

- 

0.67 

0.67 

1.13 

1999-2010 

230111 

2001BE10 

-17.126 

6.911 

-31.614 

12.758 

2.5 

19.1 

400 

0.82 

0.37 

1.79 

2001-2013* 


2005EY169 

-39.349 

17.142 

-100.220 

43.662 

2.3 

22.1 

- 

1.31 

0.23 

1.06 

2005-2014 

388189 

2006DS14 

-30.247 

13.204 

-58.690 

25.620 

2.3 

20.5 

315 

0.86 

0.34 

1.79 

2002-2014 


2004BG41 

-245.090 

110.220 

-573.400 

257.850 

2.2 

24.4 

- 

2.52 

0.61 

2.25 

2004-2011 

390522 

1996GD1 

-37.800 

17.403 

-84.884 

39.080 

2.2 

20.6 

- 

1.19 

0.35 

2.41 

1996-2014 

3361 

Orpheus 

5.220 

2.372 

12.109 

5.502 

2.2 

19.0 

348 

1.21 

0.32 

0.54 

1982-2013 

154590 

2003MA3 

-46.441 

21.733 

-96.389 

45.108 

2.1 

21.7 

86 

1.11 

0.40 

1.53 

1998-2012 

138852 

2000WN10 

15.025 

7.079 

32.239 

15.190 

2.1 

20.2 

- 

1.00 

0.30 

1.02 

2000-2013 

162361 

2000AF6 

18.534 

9.037 

33.979 

16.569 

2.1 

20.1 

- 

0.88 

0.41 

1.26 

1991-2014 

377097 

2002WQ4 

-10.360 

4.971 

-23.635 

11.342 

2.1 

19.6 

- 

1.96 

0.56 

0.92 

1950-2014 

350523 

2000EA14 

48.659 

24.108 

116.120 

57.534 

2.0 

21.1 

- 

1.12 

0.20 

2.08 

2000-2013 

7336 

Saunders 

19.567 

9.751 

53.736 

26.780 

2.0 

19.0 

600 

2.30 

0.48 

2.46 

1982-2010 

1685 

Toro 

-1.225 

0.826 

-2.732 

1.842 

1.5 

14.0 

3750 

1.37 

0.44 

1.19 

1948-2013* 

1627 

Ivar 

0.256 

0.820 

0.694 

2.221 

0.3 

12.6 

8370 

1.86 

0.40 

0.44 

1929-2014* 

433 

Eros 

-0.283 

0.361 

-0.765 

0.977 

0.8 

10.3 

23300 

1.46 

0.22 

1.40 

1893-2014* 


1 (1620) Geographos size is 5x2x1.5km. 

2 (6489) Golevka size is 350x250x250m. 


3.5. Comparison with previous results 


As already mentioned, previous studies (Farnocchia et a l.|2013bj 
Nugent et al.|20 12b) also detected a drift in semi-major axis for 
several NEAs. We compare our values with those of previous 
studies (Fig.|2|. 


Our results are in good agreement with the results previ¬ 
ously published. The difference can be mainly explained by ad¬ 
ditional observations. Fo r example, |Nugent et al. ( 2012b| ) used 
observations until 2012, [Farnocchia et al.| ( [2013b ) until 2013, 
and we used observations until May 2014. A different thresh¬ 
old for rejection may also explain some differences. For exam- 
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(101955)Bennu 
(152563)1992BF 
2005ES70 
(3908)Nyx 
(4179)Toutatis 
(2062)Aten 
(162004J1991VE 
(1862)Apollo 
(29075J1950DA 
(2340)Hathor 
(37655)lllapa 
(54509)YORP 
(1620)Geographos 
(6489)Golevka 
(283457)2001MQ3 
(2100)Ra-Shalom 
(99907)1989VA 
(230111)2001BE10 
1999MN 
(3361)Orpheus 
(154590)2003MA3 
(7336)Saunders 
(1685)Toro 
(433) Eros 


-80 -60 -40 -20 0 20 40 60 

2009BD 
2004BG41 

-600 -500 -400 -300 -200 -100 

drift (10 4 au/My) 

Fig. 2. Comparison of the drift in semi-major axis we obtained and that 
of previous studies (Farnocchia et al.[2013b| Nugent et al.|2012b|. 


i 



pie, the drift for 1999MN is sensitive to the rejection thresh¬ 
old because by using a less restrictive threshold, we have a = 
(43.2 + 10.7) x 10 -4 au/My , which is closer to the drift obtained 
by Farnocchia et al. ( |2013b| ). This may indicate a problem with 
astrometry for this object that should be remeasured. 

Finally, considering the 1-cr uncertainty, the greater discrep¬ 
ancy with previous results is for (101955)Bennu, with a differ¬ 
ence of about 4 cr. An attentive control of the astrometry and 
radar measurements do not allow us to explain this discrepancy. 

3.6. Discussion 

The drift in semi-major axis can be related to physical pa¬ 
rameters. For example, as the diurnal component is dominant, 
|Vokrouhlicky et aiTj ( |2008j i provides a relation between the drift, 
the obliquity of spin axis y, the bulk density p b , the diameter D , 
and the diurnal thermal parameter 0, i.e. 

da cosy © 

- or-r". (4) 

dt p b D 1 + © + 0.5© 2 

On the one hand, the knowledge of physical parameters can 
provide expected drift rates. For example,[Nugent et al.]( |2012a| 
used diameters and albedo from WISE data to predict drift 
rates due to the Yarkovsky effect for 540 NEAs. On the other 
hand, the knowledge of da/dt can constrain physical parameters 
D,p , or ©. In particular, if one of these parameters is already 
known, it would be helpful for the determination of other un¬ 
known parameters. 

Article number, page 6 of 1 11 1 


We also note the relation between the diameter of the NEA 
and its drift value. Diameter, absolute magnitude, and albedo of 
NEAs are available on the EARN database^] and when they are 
not available, we consider the diameter (in km) given by the re¬ 
lation with the absolute magnitude H and the albedo p v and a 
parameter K = 1329 km (Pravec & Harri sj2007| , i.e. 


DV/A = Kx 10 _H/5 . (5) 

As most of the NEAs we selected have a diameter smaller 
than 0.6 km, we can obtain statistics for the distribution of the 
diameter and the mean value of the drift (in absolute value). For 
NEAs bigger than 0.6 km, we do not have enough objects to ob¬ 
tain consistent statistics. Figure[3]represents this distribution for 
small NEAs (D < 0.6 km). The number of components of each 
bar of the histogram is also indicated. Even if the drift in semi¬ 
major axis does not depend only on the diameter (see Eq. [4]), we 
note the expected 1/D dependence of \a\. 


350 

300 





diameter (km) 


Fig. 3. Distribution of diameter and mean value of ci for NEAs smaller 
than 0.6 km. Number of component is indicated above each bar of the 
histogram and the standard deviation of the mean is represented by error 
bar. 


Another interesting result can be found in the mechanism of 
NEA transportation described in Bottke et ak| ( |2002[ l because one 
of the consequences of this mechanism is an excess of retrograde 
rotators in the NEA population. According to La Spina et al. 
p004) , 67% of the NEAs should be retrograde. 

Actually, in equation[4] the sign of d allows us to distinguish 
retrograde NEAs (negative drift) and prograde NEAs (positive 
drift). For the 43 selected NEAs from Table[4j we found 77% of 
retrogade rotators, which is close to the distribution given in La 
Spina et al j( |2004j i. 


4. The precovery context 

As a long period of observations is necessary to detect an accu¬ 
rate drift in semi-major axis, data mining and precovery obser¬ 
vations can appear very important in that context. The case of 
(152563) 1992BF studied in Vokrouhlicky et al. ( 2008) 1 is a sig¬ 
nificative example because four precoveries were available and 


http://earn.dir.de/nea/ 
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reduced with a modem catalogue to refine astrometry and then 
orbit determination. 

Nowadays, photographic plates used until the 1990s can be 
digitised and reduced with modern techniques and modern stel¬ 
lar catalogues ( [Robert et al.|2011) as the current observations on 
CCD frames. In this context, old observations can have a preci¬ 
sion equivalent to current observations. 

This section aims to investigate the case of NEAs with ex¬ 
isting precoveries for which a new reduction of these observa¬ 
tions would help to better determine the drift in semi-major axis. 
Then, by studying two specific NEAs, we highlight favourable 
periods that allow a better uncertainty of the drift, showing that 
data mining has to be considered a great opportunity in the de¬ 
tection of drift. 

4 . 1. Existing precoveries 

In this section, we identify NEAs for which precovery obser¬ 
vations have been isolated but not used for orbit determination 
because of their inaccuracy. These observations allow us to ex¬ 
tend the orbital arc for several tens of years and, in the case of 
a new reduction, they would allow for a better determination of 
the drift in semi-major axis, or at least a smaller uncertainty. De¬ 
pending on the quality of the precovery observations, we con¬ 
sider two different precisions for the precoveries with a new re¬ 
duction: 300 ma^] corresponding to accurate observations, and 
1 arcsec corresponding to low-quality observations. We deter¬ 
mine era considering both cases. In any case, the new reduction 
will decrease the uncertainty in the drift, but to select NEAs for 
which the new reduction would be helpful for the drift detection, 
we filter NEAs with the ratio a max /<Ti,, where (r u is obtained con¬ 
sidering a 300 mas precision for precovery. In fact, in the case 
in which the drift is close to its maximum a„ mx , this ratio gives 
an approximate S/N obtained with the new reduction. Table [5] 
presents 14 NEAs with a ratio a max l(Ta > 1 and indicates the 
year of the precovery, the length of additional arc (the time be¬ 
tween the first observation used in orbit determination and the 
precovery), the number of precoveries, the current uncertainty in 
the drift in semi-major axis and the expected uncertainty of the 
drift in the case the precoveries can be reduced with an precision 
of 1 arcsec and 300 mas and the ratio aWx/oy,. 

Currently, these NEAs have a bad determination of the drift, 
but a new reduction of the precoveries would help to decrease the 
uncertainty. Even if the precision of precoveries with the new re¬ 
duction is 1 arcsec (corresponding to a low-quality observation), 
we can expect a large improvement of the drift uncertainty es¬ 
pecially as the length of orbital arc is increasing. Obviously, the 
uncertainty of the drift is much smaller as we consider 300 mas 
for the precision of precoveries. 

4 . 2 . Importance of data mining 

Photographic plates used before the 1990s were used for vari¬ 
ous purposes and surveys. The field of view covered is usually 
several degrees and many objects can appear on the plates, in 
particular, NEAs leading to precovery observations. In this sec¬ 
tion, we investigate the impact of the date and the precision of 
the precovery on the drift uncertainty. To measure the impor¬ 
tance of the precovery observations, we consider the case of 
(152563) 1992BF, except with only the 1992-2011 observations 
(i.e. without the four precoveries of 1953). The idea is to under- 


7 mas is for milli-arcsec. 


stand why the 1953 precoveries are important and what happen 
if the precovery date changes. 

In this context, we consider one additional precovery obser¬ 
vation made between 1950 and 1990 and we determine the un¬ 
certainty of the detected drift with this precovery observation. 
For this, we consider various values of precision for the precov¬ 
ery observation: 3 arcsec, which corresponds to a bad astrometry 
(inaccurate stellar catalogue, bad quality of photographic plates), 
300 mas, which corresponds to mean precision of current obser¬ 
vations, and 30 mas , which could correspond to the precision 
obtained with the future Gaia stellar catalogue (see Sect. [5]). 

Figure[4]shows the uncertainty of the drift in semi-major axis 
obtained in that context for (152563) 1992BF. The grey zone 
under the curve indicates period for which <Ta is smaller than 
10 4 au/My, allowing for an accurate detection of the drift. 

If the precovery observation has an precision of 3 arcsec, 
then most of the time cru is close to 4.4 x lCD 4 au/My , which 
is the uncertainty obtained with only 1992-2011 observations. 
During short period of time (1953, 1960), the uncertainty be¬ 
comes close to2xl(T 4 au/My. These periods correspond to close 
approaches with Earth, when the geocentric distance of the aster¬ 
oid is smaller than 0.4 au. With a 300 mas precision for the pre¬ 
covery observation, we still have these favourable periods (1953, 
1960, 1966) where <r„ becomes smaller than 10 4 au/My, but it 
becomes even smaller when the precovery is farther in the past. 
This is even more the case with a precision of 30 mas for the 
precovery observation, where <Ja < 10 4 au/My , whatever the 
date of precovery observation before 1965 might be. 

Close approaches with Earth appear as good opportunities to 
improve the uncertainty of the drift in semi-major axis. Indeed, 
NEAs move faster in the celestial sphere leading to stronger con¬ 
straints on their positions and in orbit determinations. Besides, 
during close approaches, NEAs appear brighter than usual and 
can leave traces on photographic plates. 

We repeat the study for (99942) Apophis, which is a NEA 
with a very close approach to Earth on April 2029, at about 
38 000 km. Currently, we have cr* = 8.778 x 10 -4 au/My and 
very small S/N, which is obviously not significant, but provides 
an initial constraint. 

Figure [5] shows that the uncertainty of the drift could be re¬ 
duced by a precovery observation performed during favourable 
periods (1965, 1972, 1980, 1990), which corresponds to close 
approach with Earth when a geocentric distance is smaller than 
0.1 au. On these occasions, Apophis could have been imaged 
on photographic plates. Actually, we found a Beijing Observa¬ 
tory plate that covers the position of Apophis on January 1981. 
Unfortunately, because of its apparent magnitude of about 18, 
Apophis could not be detected on the plate. Nevertheless, in the 
case of a positive precovery observation at this date and with a 
precision of 300 mas, it would be possible to have a drift uncer¬ 
tainty of about 0.8 x 10 4 au/My. 

Using astrometric observations, there are two ways to de¬ 
crease the drift uncertainty: by extending the orbital arc of the 
NEA (see Fig. |T]» and during NEA close approaches with Earth. 
Figures[4]&[5]show that the drift uncertainty <r a is smaller when 
the geocentric distance is small. During these periods, NEA ap¬ 
pears brighter and during past close approaches, bright NEA 
could have been imaged on photographic plates. This is the case 
for several objects, including (152563) 1992BF, for which the 
photographic plates come from the Palomar Mountain series that 
was used for the DSS survey. This is an example of an exploited 
photographic archive, but other observatories in the world may 
have unexploited archives and could bring new observations to 
improve orbit determination and drift detection. 
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Table 5. NEAs with existing precoveries for which a new reduction will be helpful (see text). The year of the first precovery, the additional 
extended arc, the number of available precoveries, and the current and expected drift uncertainty obtained for various precision of precoveries, are 
presented. 


Asteroid 

Year of 

1 st precovery 

Additional arc 

(years) 

Number of 
precoveries 

Current ay 

Expected oy 
(1 ") 
(10~ 4 au/My) 

Expected oy 
(300 mas) 

&max ! O' ci 

162421 

2000ET70 

1977 

23.050 

2 

16.8 

1.7 

0.7 

8.04 

4179 

Toutatis 

1934 

42.292 

2 

0.4 

0.4 

0.3 

4.43 

138175 

2000EE104 

1999 

0.923 

5 

11.6 

5.8 

4.5 

3.46 

7350 

1993VA 

1963 

22.563 

3 

7.7 

2.0 

1.5 

2.37 

66400 

1999LT7 

1987 

12.134 

2 

19.1 

7.0 

4.5 

2.29 

5797 

Bivoj 

1953 

26.253 

1 

4.3 

3.6 

3.4 

2.28 

350751 

2002AW 

1991 

10.776 

2 

72.9 

24.2 

9.8 

2.01 


1999FA 

1978 

19.855 

1 

127.7 

10.3 

9.7 

1.96 

87024 

2000JS66 

1979 

21.328 

2 

32.4 

11.0 

4.8 

1.54 


2004SV55 

1983 

20.997 

5 

576.5 

4.7 

3.9 

1.31 

7341 

1991VK 

1981 

9.886 

1 

3.3 

3.2 

2.4 

1.27 

1943 

Anteros 

1968 

4.767 

2 

1.9 

1.9 

1.7 

1.05 

5604 

1992FE 

1976 

8.983 

2 

10.0 

5.5 

3.5 

1.01 

4450 

Pan 

1937 

50.631 

2 

10.8 

4.7 

3.7 

1.01 





1950 1955 1960 1965 1970 1975 1980 1985 1990 

date 


Fig. 4. Uncertainty of the drift in semi-major axis related to the date 
of the precovery observations for (152563) 1992BF, with a precision of 
3 arcsec, 300 mas , and 30 mas from top to bottom, for the precovery 
observation. Unit is 10 _4 au/My for cr b in ordinate. 





Fig. 5. Uncertainty of the drift in semi-major axis related to the date 
of the precovery observations for (99942) Apophis, with a precision of 
3 arcsec, 300 mas, and 30 mas from top to bottom, for the precovery 
observation. Unit is 10 _4 au/My for oy in ordinate. 


5. The Gaia mission context 


The space mission Gaia was launched successfully on December 
2013. The aim is to have a 3D-map of our Galaxy. In the case of 
solar system objects, the satellite will enable us to map thousands 
of main belt asteroids and NEAs down to magnitude 20. Gaia 


will observe several times each NEA for five years of the mis¬ 
sion. [Mouret & Mignard ( 2011 [ ) identified 64 promising NEAs 
for which Gaia observations could help the detection of a drift in 
semi-major axis with an uncertainty better than 5 x 10 4 au/My. 
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Beyond the Gaia observations themselves, one of the best ex¬ 
pected improvements will be the Gaia stellar catalogue allowing 
for orbital improvement and detection of drift in asteroid semi¬ 
major axis. The Gaia catalogue ( |Mignard et al.||2007] ) will pro¬ 
vide unbiased positions of a billion stars down to magnitude 20 
and with a precision depending on magnitude (7 yuas at <10 mag; 
12-25 /ras at 15 mag and 100-300 /tas at 20 mag). 

With such a precision, new processes of reduction will be 
necessary, such as propagation models to third order for star 
proper motions, differential aberration, atmospheric absorption, 
color of stars, etc. Finally, by including these corrections and 
considering the limitation due to the seeing, an uncertainty of 
10-20 mas can be expected on the position of any asteroid re¬ 
duced with the Gaia stellar catalogue. 

Moreover, the Gaia stellar catalogue will enable us to reduce 
old photographic plates. The reduced position will be affected 
by other errors or limitations, including the distortion or instru¬ 
ment deformations; the granularity of the photographic plate; the 
digitisation; the exposure time; and the proper motions of stars. 
Depending on the magnitude and the epoch, the error would be, 
for example, about 12 mas for 20-mag star and 1.5mas for a 16- 
mag star at epoch 1950 ( jMignard|2013[ ). 

Another problem could be raised from a lack of meta-data in 
time observations. The time of the observations could be given 
with a low precision (to the minute, for example), or an error 
clock in the time of observation could possibly lead to inaccurate 
or even worse unusable observations. 

In this study, we do not consider clock error, but only errors 
due to digitisation, photographic plates, or instruments. Assum¬ 
ing a good plate quality, fine grain emulsion, and optimal plate 
measuring, [Zacharias| ( j2013"| estimated an error budget of about 
10-20 mas for an astrometric position extracted from digitised 
plate. For these reasons, we assume that the astrometric preci¬ 
sion of the new reduction with Gaia stellar catalogue will be 10 
mas for observations since 1990 (for mainly CCD observations), 
30 mas for 1950-1990 observations, 50 mas for 1920-1950 ob¬ 
servations, and 100 mas for observations before 1920 (mainly 
photographic plates). 

To measure the improvement of the drift detection with the 
Gaia stellar catalogue, we choose several assumptions as to new 
reduction methods that make sue of the Gaia catalogue. It seems 
obvious that all of the observations could not be reduced again, 
thus we assume that only a few would be reduced. As mentioned 
in Desmars et al. (20131, it appears that in such a case, the best 
solution for orbit determination is to reduce the first and last ob¬ 
servations rather than a random part of them. Actually, the first 
and last observations, by virtually increasing the length of orbital 
arc, bring more important constraints on the orbital motion of the 
asteroid. In that context, we consider five cases (also considered 
in |Desmars et al.|2013| l : 


- Case 1: we examine current observations with their current 
uncertainty; 

- Case 2: we assume that the first and the last observations can 
be reduced with the Gaia stellar catalogue; 

- Case 3: we assume that the first five and last five observations 
can be reduced with the Gaia stellar catalogue; 

- Case 4: we assume that the first ten and the last ten observa¬ 
tions can be reduced with the Gaia stellar catalogue; 

- Case 5: we assume that all observations can be reduced with 
the Gaia stellar catalogue. 


Table [6] presents the distribution of drift uncertainty for the 
1629 numbered NEAs according to various cases of new reduc¬ 
tion with Gaia catalogue. 


By considering different thresholds for the uncertainty, we 
note that the Gaia catalogue will be helpful to get accurate de¬ 
tection of drift in semi-major axis. Even in the minimal case 
where only two observations can be reduced with the Gaia 
catalogue (the first and last observations), 31 NEAs will have 
iTa < 10 4 au/Mv, whereas only eight NEAs currently reach this 
level of uncertainty. The more we reduce observations with Gaia 
catalogue, the more the quality of the drift detection will im¬ 
prove. In the extreme case where all observations could be re¬ 
duced with Gaia catalogue, 688 (more than 40% of numbered 
NEAs) will have cru < 10~ 4 au/My. 

Considering the maximum value of the drift, the ratio 
a maxima also helps to measure the impact of the Gaia catalogue 
in the detection of drift. Table [ 7 ] presents the distribution of ra¬ 
tio for the different cases of reduction. As mentioned in previous 
section, the ratio can provide an approximation of the expected 
S/N in the case of the drift is close to its maximal value. So 
Umax/va > 3 could lead to reliable detection of the drift. Cur¬ 
rently, only 15 objects reach this ratio, but with the minimal 
case with only two observations reduced with Gaia catalogue, 
70 NEAs will have ci max lcra > 3. Reducing 161 and 223 NEAs 
with 10 and 20 observations, respectively, using the Gaia cata¬ 
logue, will fulfill the same condition. Finally, in the extreme case 
where all observations could be reduced again, 938 NEAs (more 
than 57%) will have a ratio greater than 3. 

The Gaia catalogue will be very useful for orbit determina¬ 
tion and for accurate detection of drift in semi-major axis, even 
with a few number of new reduction with Gaia catalogue. Of 
course, when the Gaia catalogue will be available in 2022 (or 
even with intermediate releases during the mission), new obser¬ 
vations will be done and the drift uncertainty will be even better 
for most of NEAs. 


6. Conclusion 

The Yarkovsky effect produces a drift in the semi-major axis of 
NEAs. With a precise process of orbit determination, we have 
detected consistent drift for 46 NEAs. The drift value we mea¬ 
sured are in agreement with values obtained by statistical meth¬ 
ods. Moreover, the robustness has been tested for three represen¬ 
tative NEAs. The drift in semi-major axis can be related to other 
physical parameters, such as the diameter, the obliquity, or the 
thermal parameter. The knowledge of the drift allows us to con¬ 
strain these physical characteristics. For example, the prograde 
or retrograde rotation of the NEA can be determined by the sign 
of a. The ratio of 78% of retrogade rotators is consistent with the 
expected ratio explaining the mechanism of NEA transportation. 

We have also investigated the importance of precovery ob¬ 
servations in the accurate detection of drift. In particular, we 
have shown that the drift uncertainty can be deeply improved by 
reducing only one precovery with modern techniques and mod¬ 
ern stellar catalogue. The improvement is more important during 
close approach with Earth, which usually corresponds to a good 
period of visibility for the NEA (smaller apparent magnitude). 
During previous close approaches, some NEAs may have been 
imaged on photographic plates, and thus could provide accurate 
observations in order to refine orbit determination or, even more 
importantly, to help to detect a drift. We have already identified 
14 NEAs with unused precoveries because of bad astrometry that 
can help to refine the detection of the drift. 

Finally, we have highlighted and quantified the impact of the 
Gaia catalogue in the drift determination. With this catalogue, 
future and also past observations could be reduced with very 
high precision (a few tens of mas). If we reduce a few current 
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Table 6. Number and percentage of numbered NEAs for which the uncertainty of the drift will be smaller than representative values with different 
assumptions of new reduction with Gaia catalogue. 


Case 

Number of observations 
reduced with Gaia catalogue 

CTti < 10 

3 au/My 

CTa < 10 

4 au/My 

(Ta < 10 S au/My 

1 

0 

107 

(6.6%) 

8 

(0.5%) 

1 (0.1%) 

2 

2 

346 

(21.2%) 

31 

(1.9%) 

1 (0.1%) 

3 

10 

494 

(30.3%) 

93 

(5.7%) 

4 (0.2%) 

4 

20 

644 

(39.5%) 

149 

(9.1%) 

6 (0.4%) 

5 

all 

1442 

(88.5%) 

688 

(42.2%) 

91 (5.6%) 


Table 7. Distibution of ratio a max la il with different assumptions of new reduction with Gaia catalogue. 


Case 

Number of observations 
reduced with Gaia catalogue 

[0,1] 

[1,3] 

&max /^a 

[3,10] 

[10,20] 

>20 

1 

0 

1567 

47 

14 

0 

1 

2 

2 

1423 

136 

61 

7 

2 

3 

10 

1291 

177 

117 

29 

15 

4 

20 

1152 

254 

135 

56 

32 

5 

all 

319 

372 

461 

196 

281 


observations with the Gaia catalogue, we can deeply increase the 
uncertainty in the drift detection. For example, with only two ob¬ 
servations reduced again with the Gaia precision, 31 NEAs will 
have (Ta < 10 4 au/My whereas only 8 NEAs currently reach this 
uncertainty. The Gaia catalogue opens up new perspectives in 
the drift detection. 

The reduction of old photographic plates of solar system 
objects in the context of the Gaia mission is the topic of the 
NAROC^] project (Arlot 2013 ). This project investigates natural 
satellites and asteroids (Transneptunian objects, NEA, or Main 
Belt asteroids) in order to refine the dynamics of these objects. 
The detection of the drift in semi-major axis of NEAs is an im¬ 
portant part of this project. 
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Appendix A: Maximum drift 

To highlight unreliable detection, we introduce an estimation of 
a maximal drift in semi-major axis related to the absolute mag¬ 
nitude H. |Vokrouhlicky| ( | 1 998[ > then |Vokrouhlicky et al.] ( |2Q0Q] > 
gave an relation between the drift in semi-major axis due to di¬ 
urnal component of the Yarkovsky effect and several physical 
parameters, valid for circular orbit. We find 



8cr 

9/7 


®(a) 


Gsind 

ITT 


cos y 


(A.l) 


where a is the absortivity (complement of Bond albedo A), 
n is the mean motion, y is the obliquity, and <t>(a) is the stan¬ 
dard radiation function (inversely proportional to the diameter 
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D, the bulk density p, and t he semi-maior axis a). The expres¬ 
sion is fully detailled in Vokrouhlicky et al. ( 2000|> and can 
be approximated by 


-0.5© 


G sind 


1 + A 1 + © + 0.5© 2 


(A.2) 


where © is the diurnal thermal parameter. 

By using p = ira 3 and numerical values for <I> function, we 
have 


( 



6.412 x 10 3 


a 0.5© 

Dp yfa 1 + © + 0.5© 2 


cos y, 


(A.3) 


with the convenient units, a in au and ^ in au/My. The 
expression ^ is maximum when y = 0 or 180 degrees and 

l+gruTse 2 reaches its maximum value 0.207 for © = V2. 

Finally, 


= 1.327 x 10 


Dp yja 


(A.4) 


To obtain a simple relation, we can consider a — 1 au, a mean 
density p = 2600 kg.m 3 et a mean albedo p v = 0.154 for NEAs 
( |Chesley et al.|2002). The Bond albedo will be A — 0.060 using 
the relation from|Bowell et al.|(|1989|) A = (0.290 + 0.684G)p,. 
(with G = 0.15) and a = 0.939. 

The diameter can be approximated by the relation (Pravec &| 
|Harris|2007| l 


D \[jp, = 1329km x 10^ /5 


with H the absolute magnitude. 

Finally, the maximum drift can be approximated with the ab¬ 
solute magnitude H as 

= l0 a2H - 6 - 85 au/My. (A.5) 

max 
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